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Abstract 

We extend a recent numerical calculation of the pion electromagnetic form factor F n (Q 2 ) in 
holographic QCD to study two important issues regarding the behavior of fields in the bulk. First, 
we show that using a chiral symmetry-breaking field formally satisfying the boundary conditions of 
the "soft-wall" model changes numerical results very little from the earlier calculation that ignores 
these constraints. Second, we use a background field that interpolates between "hard-wall" and 
"soft-wall" models to obtain an improved model that reproduces the desirable phenomenological 
features of both. In all cases, F w for large Q 2 is shallower than data, an effect that can be cured 
by relaxing the fit to one of the static observables, particularly the decay constant f n . 
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I. INTRODUCTION 



One of the most interesting developments in QCD in recent years has been the applica- 
tion to hadronic physics of the gauge/gravity correspondence jl] between strongly-coupled 
gauge theories and weakly-coupled gravity on curved spacetime backgrounds. This corre- 
spondence is most firmly established between Af = 4 supersymmetric Yang-Mills theories 
(which are conformal field theories [CFT]) and a 5-dimensional anti-de Sitter (AdS) grav- 
ity background — hence the moniker AdS/CFT. However, one may conjecture that a similar 
connection holds for any strongly-coupled gauge theory possessing an approximate conformal 
symmetry, a property that holds for QCD in its high-energy limit. The terms "AdS/QCD" 
or "holographic QCD" refer to studies based upon the premise that QCD belongs to this 
class and possesses a suitable gravity dual in 5D space. Properties such as confinement 
and chiral symmetry breaking arise in many models built in this way; indeed, in the "top- 
down" approach one starts with a string theory and chooses a gravity background that 
reproduces such basic QCD features. Much more amenable to phenomenological analysis is 
the "bottom-up" approach, in which one begins with the QCD observables and determines 
what gravity backgrounds reproduce them most successfully. 

The fields of our 4D universe occupy only one surface of the 5D AdS space, whose metric 
is given by 



ds 2 = g MN dx M dx N = — {ji^dx^ dx v — dz 2 ) . (1) 
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Here the full nontrivial 5D metric g MN obtained from Eq. ([1]) is distinguished from the 4D 
Minkowski metric r]^ = diag(+, — , — , — ). The ultraviolet (UV) limit of QCD is represented 
by fields living near the AdS singularity z = e — > (the "UV brane"), suggesting the as- 
sociation of the extra "bulk" coordinate z with inverse momentum scales: Q ~ 1/z. The 
gauge/gravity correspondence states that every CFT operator 0(x) is associated with a bulk 
field *&(x, z) of given quantum numbers uniquely determined by its value e) on the UV 
brane, which explains the origin of the usage "holographic". Additionally, the global QCD 
symmetry of isospin associated with the two light quark flavors is promoted to a gauged 
SU(2) symmetry respected by the bulk fields. 

From the point of view of hadronic physics, the original J\f=4 super- Yang Mills theory 
is inadequate due to being exactly conformal and therefore lacking asymptotic S-matrix 
particle states. In AdS/QCD one breaks the conformal symmetry (and introduces a mass 
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scale) by impeding the ability of the fields *$?(x, z) to penetrate deeply into the bulk, leading 
to an explanation of confinement and more generally constraining the model's infrared (IR) 
behavior. The most straightforward realization uses a hard cutoff at a particular value 
z = zq ~ 1/Aqcd (the "IR brane") at which point appropriate boundary conditions are 
imposed upon ^Sf(x,z), thus defining the so-called "hard- wall" model |2j. However, despite 
its economy, the hard- wall model leads to results incompatible with the expectations of linear 
confinement: For any model in which \[f(x,z) penetrates to a limited depth into the bulk, 
the eigenmodes of ^ (with mass eigenvalues m n ) represent hadronic states all carrying the 
same quantum numbers, producing the AdS/QCD version of Regge trajectories. However, 
QCD with linear confinement is expected to follow the trajectory m 2 ~n, while hard- wall 
models predict rn? n ~n 2 . To repair this shortcoming requires the introduction of a "soft-wall" 
gravity background with an exponential decrease ~ e~ K2z2 in the action for large z |4|, in 
which case k serves the role of Aq CD . 

Studies of hadronic properties in AdS/QCD have been very popular in recent years; 
to name but a few that share the same spirit as this work are examinations of hadronic 



spectra ^Jfi, fl EJs], I^L lj 
breaking (a, U, U, 



11| , the couplings of hadrons in the presence of chiral s ym metry 



hadronic form factors 
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15 L the quark- quark potential [161 ]. 4-quark operators 17], and 
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23(| . In this paper we continue a study of the 
pion electromagnetic form factor F n initiated by the present authors in Ref. 24J], using the 
treatment of chiral symmetry-breaking effects (parametrized by a light quark mass m q and 



chiral condensate a) explicitly incorporated into the Lagrangian, in the manner deve 



by the model of Ref. [5j. F n has also been considered in two other recent papers 22 



oped 



23|; 



uses a model that does not incorporate chiral symmetry breaking explicitly, while 



Ref. b2] 

Ref. [23] also uses the formalism of Ref. |5| and has a considerable overlap [37j with the 



hard- wall calculations of Ref. 



24| ; indeed, it also proves a number of analytical results that 



hold in the limit m q = 0. However, Ref. 23J] does not present calculations in the soft-wall 
model due to objections that we address below. 



While results for most hadronic quantities tend to agree surprisingly well with the resu 



m 



ts 



of the hard-wall model (and generally somewhat less well for the soft-wall model 
observables in the vector meson sector tend to depend somewhat less sensitively on the 
precise nature of the IR boundary condition. However, the axial-vector sector and especially 
the pion depend upon chiral symmetry-breaking in a much more direct fashion 4j ; exploring 
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this sensitivity was an original motivation of the studies in Refs. 23|, \24 \. 

Our previous work [24| numerically calculated F n (Q 2 ) for spacelike values Q 2 > of 
momentum transfer, as well as couplings fiUn)^ obtained from the timelike region, in both the 
hard-wall and soft-wall models. Our results showed that both models predict Q 2 dependence 
for F n (Q 2 ) that tends to be too shallow for both models (in particular, predicting a value 
too small for the charge radius (r 2 )), but significantly worse for the soft-wall case. Better 
agreement with F n (Q 2 ) data could be achieved only by loosening the fit to other observables, 
especially by decreasing the decay constant /„-. 

However, a small swindle is introduced in applying the methods of Ref. j^] to the soft- 
wall model 4j: The background scalar field X{z) whose vacuum expectation value X (z) 
provides m q and a via its small- z asymptotic value (2X Q (z) = m q z + az 3 ) satisfies in general 
a 2nd-order differential equation whose solutions are Kummer functions, only one of which 
[zU(^, 0, k 2 z 2 ), to be precise] satisfies the necessary finiteness boundary condition of the 
soft-wall model as z — > oo. But this particular function carries a unique proportionality 
between the z and z 3 terms of its Taylor series, implying an unphysical fixed ratio between 
m q and a. Such an obstacle is easily overcome by including higher-order terms in X(z) in 
the potential that do not greatly modify the small- z behavior; however, from the operational 
point of view, in Ref. 24( we simply observed that such large- z modifications are numerically 
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heavily suppressed by the background factor e~ K z . One of the goals of the current work is 
to present a calculation using a field Xq(z) that satisfies both constraints, finite for large z 
and allowing independent m q and a values. Numerous functional forms for the field Xq(z) 
satisfy this requirement, but for sake of definiteness we choose 



2X (z) = (m q z + az 3 ) 



1 — exp 



A. 



K±Z 4 



B c exp 



4 K 2 z 2 



(2) 



We discuss the motivation for this particular choice in Sec. IHIl but point out that numerical 
simulations for several other similar forms do not significantly alter our conclusion: Using 



a functional form satisfying a 
the naive soft-wall fit of Ref. 



1 required asymptotic behaviors does not significantly change 



24]. 



The other goal of this paper is to develop a model that combines the numerical successes 
of the hard-wall model with the Regge trajectories predicted by the soft-wall model. In 
the closing statements of Ref. [24j we presented a background field with precisely these 
properties. Inspired by the old Saxon- Woods generalization of the hard-cutoff nuclear density 
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model, we proposed 

which has a drop-off at z = zq (and equals | there), but decreases as e~ x2z2 for large z. 

Our results indicate that neither of these modifications substantially alter the outcome of 
Ref. [24|. The soft-wall model test function X (z) of Eq. (El), although formally possessing the 
correct asymptotic behaviors, gives an optimal fit to F n (Q 2 ) no better than that in Ref. 24]. 
Similarly, as expected, the adoption of Eq. (j3J) for the background field interpolates results 
between hard- and soft-wall models. This latter result may be somewhat surprising since 
Eq. ([3]) contains two independent mass parameters (1/zq and A) rather than the one of the 
strict hard- or soft-wall models. Since our previous work 24| showed that both hard- and 
soft-wall model predictions for F n (Q 2 ) can be improved significantly by loosening the fit 
to just one observable, one might expect that Eq. ([3]) with its extra fit parameter allows 
sufficient latitude not only to give both hard- and soft-wall behaviors, but also to improve 
the global numerical fit to F n (Q 2 ) and other low-energy observables. As seen below, the 
former goal is achieved but the latter is not. Taking these results together, it appears 
that achieving a better simultaneous fit to F n (Q 2 ), f n , m n , m p , f p , and other low-energy 
observables requires a more general treatment of chiral symmetry breaking than the one 
developed in Ref. [5fl. 

This paper is organized as follows: In Sec. [Til we recount the formalism of Ref. [5] relevant 
to our calculations and give expressions for F n (Q 2 ) and couplings in terms of AdS/QCD 
mode wave functions. Section [TTTI presents the results of numerical simulations of F V (Q 2 ) and 
other low-energy hadronic observables, and Sec. IIVI summarizes our results and concludes. 
Details of the numerical procedure are discussed in the Appendix. 



II. FORMALISM 



The content of this section is nearly identical to that in Ref. 24j, but is included for 
completeness. The full 5D action 5j used in this work reads 



S = Jd\ C -*« y/g Tr j \DX\ 2 + 3|X| 2 - -^(F 2 + Fr)\ 



(4) 



where g = \ detg MN \ is obtained from the metric in Eq. ([T]), and e ®^ represents a back- 
ground dilaton coupling. The original holographic QCD calculation in Ref. f| uses the 



hard-wall (step function) background e 



-*(*) 



H(z — z), while Ref. |4[ defines the soft-wall 



model by e 



', which as noted above, reproduces the traditional Regge trajectory 



behavior for the mesons. In this paper we also consider the Saxon- Woods form for e^®^ 
defined in Eq. ([3]). The scalar field X(z) [actually (2/ z)X(z)] is the holographic partner of 
the quark condensate q R q L , and its vacuum expectation value X (z) = \v(z) = ^m q z + ^az 3 
gives rise to explicit and spontaneous chiral symmetry breaking. This expression for X (z) 
is exact in the hard-wall model, while it holds only for small z in the soft-wall model. 

The chiral gauge fields A a L R (holographic partners to the bilinears q L jff^^qj^ R ) enter 
through D M X = d M X-iAfX + iXAf, A% R = Af^ R t a , and F%g = d M A 1 [ R ~d N A^ R - 
i[A% R ,A% R }. The polar V and axial A gauge fields are V M = \(A% + A%) and A M = 
\{A%-A%), in terms of which D M X = d M X-i[V M , X}-i{A M } X}, F™ N = d M V N -d N V u - 
i ([V M , V N ) + [A M , A N ]) , FX IN = d M A N -d N A M -i ({V M , A N ] + [A M , V N ]), and 



S = |^e-^v^Tr||DX| 2 + 3|X| 2 -^(^ + Fi)} . 



(5) 



We use the axial-like gauge V z (x,z) = 0, A z (x, z) = 0, so that the associated sources may 
be expressed as divergences over just the usual four spacetime dimensions, d^V^ = (since 
isospin is conserved) and d^A^. A^ can be further decomposed into a transverse (divergence- 
less) piece v4 M _i_ and a longitudinal piece <p: A^ = A^± + d^ip. The pion field ir a appears 
through X = X exp(2iir a t a )] n a is dimensionless and related to the canonically-normalized 
pion field n a of chiral Lagrangians via ir a = 7c a / f n , with ^ = 93 MeV. 

The equations of motion obtained from Eq. ([5]) for the fields ^(q, z) {i.e., all except for 
X are Fourier transformed with respect to the 4D coordinates x) read 



a 



o,v;: \ +- 







z 



2 p -*( Z ) 

0,X ) + - Al 



9, 



-*(*) 



z 



3s X 



z 



- q 2 d z ^ + 9 -t^l dz ^ 



-ax r 



+ 



(6) 

(7) 

(8) 

(9) 
(10) 



6 



The gauge/gravity correspondence itself provides the first prediction, the value of the 5D 



gauge coupling g 5 Q]. Substituting Eq. 
boundary term for the action: 

1 



into the F v portion of Eq. (j3J) leaves only the 



S 



2g. 



d x 



-*(*) 



(11) 



The vector field V?(q,z) can be resolved as V£(q, z) = V(q, z)V"(q), where V®(q) is the 
Fourier transform of the vector current source = q^pt a q at the UV boundary z = e, 
and V(q,z) (the "bulk-to-boundary propagator") is normalized to V(q,e) = 1. Due to the 
isospin conservation constraint q^V 1 = 0, one may replace V®V^ a with V^V^W LV '5 ab and 
IP 1 * = r]^ v — q^q v jq 2 i and then the usual quadratic variation of the action with respect to the 
source V produces the vector current two-point function: 



d 4 xe^(j;(x)J>)> = 5 ab n^S y (g 2 ), 

e -*Wd z V(q,z) 



gi 



(12) 
(13) 



from which one finds, matching to the QCD result for currents J M normalized [3] according 
to the prescription of Q], 

2 I 271 "' A 2 / n „X 

^5 = ^-^4vr 2 . (14) 

An analogous calculation in the axial sector relates the bulk-to-boundary propagator 
A(q, z) to the 7r decay constant /„-: 



2 ___1 d z A{0,z) 

J n 



9l 



z 



(15) 



from which one sees that the relative choices of normalization for /„. and g$ are correlated. 

The normalizable eigenstates of Eqs. (EI) (J9j) correspond to towers of hadrons of the same 
quantum numbers as the parent fields. Since large N c is intrinsic to this procedure, the 
mesons have narrow widths and the spectral decompositions of self-energy functions such as 
Sy are sums over poles: 

f 2 

J n 



M? 2 ) = £ 



n=0 



q 2 -Ml" 



(16) 



where M n are the mass eigenvalues and /„ are the decay constants of vector modes ip n {z) 
normalized according to 



dz- 



-*(*) 



1pm{z)i) n (z) = 8 n 



(17) 



In order to obtain the pion form factor F n , one expands the action Eq. (jSJ) out to cubic 
order in the fields. Since the pion field is related [Eq. Qj] to the longitudinal mode d^ip 
of A* 1 , one must identify not only Virn terms, but also VAA and VAtt. A straightforward 
calculation (that uses the equation of motion for V to simplify the result) produces the Vw 
terms: 



^AdS — ^abc I d X I dz t ^ ^ 



95 z 



+ V -!^L (#V> _ Q^a^ yb _ ^ 



(18) 



where the z integration range is [0, oo). Reference jfj uses this action (but only with e - *^ = 
1) to obtain the Vtctt couplings [Eq. (125]) belo w], with the caveat that terms cubic in Fy t A 
have not been included. In fact, we have shown [24J that no such terms contribute to the Vw 
coupling. The 3-point correlator is obtained from a straightforward variation of Eq. (JTHl) : 

(J:( Pl )J$ b (q)JZ(-p 2 )) =e abc F(plplqi)(p 1 +p 2 r i (27rr5^(p 1 -p 2 + q) , (19) 

where, again recalling the narrowness of resonances, one may express the dynamical factor 
F(pf,p 2 , q 2 ) in terms of transition form factors: 

and F n k(q 2 ) correspond to form factors for n^k transitions. The pion form factor F n (q 2 ) is 
then defined as the ground-to-ground pseudoscalar meson transition to the vector current: 

FAd 2 ) = F u (q 2 ) = [ dz e"»W { ^[^(z)] 2 + ^ Hz) - ^)] 2 } , (21) 

J Jn K.95 Z Z J 

whose normalization F n (0) = l [using V(0, z) = l] is guaranteed by the canonical normaliza- 
tion of the pion kinetic energy term in the 4D Lagrangian. The origin of the z dependence of 
Eq. ([2T]) is recognizable from Eq. ( |T8i) . The pion is the ground-state solution of Eqs. (jEJ)-© 
for ir a subject to appropriate boundary constraints: Neumann at z = Zq for the hard-wall 
model, finite as z^oo for the soft-wall model. The solution to Eq. ([6]) for source V in the 
soft-wall model [<&(z) = n 2 z 2 ] for spacelike momentum transfers q 2 = —Q 2 < 19j, |22j, l24j . 



used in our first calculation, includes the Kummer (confluent hypergeometric) function U: 

v( q , z) = r(i + q 2 /ak 2 ) u[q 2 /ak 2 , o, M 2 ] . (22) 



Note that for large z, Eq. (122]) falls as (z 2 )~® 2 / 4k2 . For our second calculation [e - *^ as in 
Eq. (jBD], z) must be obtained numerically from Eq. ([6]). In both cases the solutions 
satisfy the boundary conditions V(q, e) = 1, V(0, z) = 1. In the timelike region, V(q, z) may 
be expanded as 

n=l ^ n 

When numerical solutions are required, ip n (z) and / n are obtained by solving Eq. at its 
poles and properly normalizing using Eq. ffTTj) . In either case, substituting Eq. (123]) into 
Eq. (121 p gives the timelike pion form factor as a sum over vector meson poles: 



fn9 

n=l 

where g n7Tn is given by 



^« 2 > = -£#riJi- < 24 > 

n=l y « 



<W = § fdzMz)e- Hz) \^ldMz)] 2 + ^k(z)-rtz)f) ■ (25) 

Together, Eqs. (I2T1) and Eqs. (j2l|) - (|25p . with V and ^ n obtained from either Eq. (1221) or the 
numerical solution to Eq. (jHJ), provide a complete expression for F n (q 2 ) at all values of q 2 . 



III. RESULTS 

For our first result, we consider the consequences of modifying the solution for the vev of 
the bulk field X(z), as stated in the Introduction and repeated here for convenience: 



2X0(2;) =v(z) = (m q z + a z 3 ) 



exp 



K 4 Z A 



+ ^ex P [- J J- 5 ]. (2) 



This particular expression is chosen because of the interesting large- z asymptotic form of 
solutions of Eq. ( flOl) for Xq(z) [J]: In the hard- wall case the independent solutions are 
precisely z 1 and z 3 (so that the asymptotic solution is also the full solution in that case), 
while for the soft-wall case the exact solutions are the Kummer functions z 3 M(|, 2, k, 2 z 2 ) 
and zU(l,0, k 2 z 2 ). Of the latter pair only the second is finite as 2:^00, and indeed behaves 
as exp(— 3/4k 2 z 2 ). As indicated in the Introduction, Xq must have corrections arising from 
higher-order terms in the potential to allow independent m q and a parameters at low z, which 
in turn would appear as corrections to Eq. ( 1101) . This argument originated with the original 



soft-wall paper [4], and since a full explanation that motivates our choice Eq. (j2J) requires 



understanding the original argument in detail, we repeat it here. Inasmuch as e - *^ vanishes 
as z — ► oo in the soft-wall model, the equation of motion Eq. (flOl) for X is guaranteed to 
possess a solution that tends to a constant as z^oo. Since the e - *^- 1 factor also multiplies 
any higher-order potential terms in X that might be introduced into the action Eq. (jSJ) 
and hence into Eq. (fTOI) . the full equation of motion continues to possess a solution that 
approaches a constant as 2-> oo. Moreover, for such a solution with Xo(oo) = constant, 
all multiple-derivative terms appearing in Eq. (1101) and its generalizations become negligible 
compared to the lone single-derivative term contained in the first term of Eq. ( FlOl) . The 
constant X (oo) therefore acts as the sole integration constant in this effectively lst-order 
differential equation, and thus determines the (in general nonlinear) relationship between its 
solution's z 1 and z 3 coefficients, namely, the ratio of m q to a. By adjusting the higher-order 
terms in the potential, one changes Xo(oo), thus altering the numerical ratio between its 
coefficients near z = 0, thereby allowing one to obtain the desired ratio of m q to a. 

With these features in mind, let us return to the motivation for the choice Eq. (J2J). The 
asymptotic form exp(— 3/4k 2 z 2 ) already appears for the explicit solution to Eq. (flOl) and 
therefore is chosen for the z—>oo limit of Eq. (j2J). The functional forms exp(— A/z n ) with 
A,n>0 are especially interesting because they not only rapidly approach 1 as z— >oo, but 
also sluggishly depart from as z is increased from 0, due to their nonanalyticity about 
that point. Thus the small-z asymptotic form of Eq. (j2J), m q z+az 3 , is precisely as required. 
These properties make Eq. an excellent choice for a test function satisfying the desired 
asymptotic behaviors. Note that we do not exhibit a specific potential with this solution to 
its equation of motion, since Eq. (j2J) satisfies all the requirements placed upon it through 
its solution. Moreover, Eq. (j2j) gives results for F n (Q 2 ) and other observables in particular 
ranges of its parameters A c and B c that agree numerically with the original soft-wall model 
(which, recall, is the exact solution for a particular potential), thereby providing additional 
(a posteriori) support for this choice. In particular, the precise form of X Q at intermediate 
values of z — ignored in the construction of Eq. (JSJ) — is apparently not important in fitting 
observables. 

As mentioned above, the numerical change from the results of Ref. 24j (the original soft- 
wall model curve in Fig. CP due to using Eq. (j2J) turns out to be quite minimal in the best-fit 
case. For sake of definiteness, let us choose B c — 1 and A c — 1. We find these values give a 
best fit to F n (Q 2 ) data quite close to (but slightly worse than) that of the original soft-wall 
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TABLE I: Comparison of soft-wall model to modified e *( z ) with Xzq = 1; values in MeV (except 



for g p7T7V ). 

Observable Experiment Soft-wall Xzq = 1 

m n 139.6±0.0004 [25] 139.6 139.6 

m p 775.5±0.4 [25] 777.4 779.2 

m ai 1230±40 [25] 1601 1596 

U 92.4±0.35 [25] 87.0 92.0 

/ p 1/2 346.2±1.4 [26] 261 283 

/a 1 / 2 433±13 [27, 28J 558 576 

g p7T7T 6.03±0.07 [25] 3.33 3.49 



model obtained in Ref. [24J — an encouraging start. Keeping A c — 1 and varying B c , we find 
that values of B c < 1 (including I? c <0) give an F n prediction that moves even further away 
from the data (i.e., is even more shallow) than the original soft-wall prediction. On the 
other hand, for B C >1 the prediction moves closer to F n data, but still is only as good as the 
original soft- wall prediction at B c « 2. For B c >2 the prediction again moves away from the 
F n data. Next, if we keep B c = l and vary A c , we find that for < A c < 1 the prediction moves 
away from the F n data, but for A c > 1 the prediction improves only marginally. Interestingly, 
letting both B c and A c vary always appears to permit a set of parameters that fit the data [F n 
and the static observables m n , m p , m ai , f n , f p ', fl[ 2 , and g p7T7T ] reasonably well, but only if 
B c <2. However, these fits are never superior to the original soft-wall fit. We conclude that 
the original fit in Ref. 2^], although formally violating the AdS/QCD boundary conditions, 
is actually quite stable when the proper boundary conditions are imposed. 

That the new fits never improve upon the original model is quite easy to explain nu- 
merically. First note that V(q, z) in Eq. (122]) drops from 1 at z — to smaller values for 
larger z with increasing rapidity as Q 2 increases, a feature not affected by our modification 
to Xq(z). The effect of the new Xq(z) is reflected only in the rest of integrand in Eq. ( !2~TI) . 
whose integral [i.e., for V(q, z) — > 1] is normalized to 1. One expects that the integrand 
favors larger z in the original soft-wall case because of the direct z 3 dependence of Xq(z) 
than in the modified case, since X (z) — > B c as z—>oo. Since V(q,z) decreases with z for 
all Q 2 >0, F n (Q 2 ) is also lower in the original model for all values of Q 2 >0. Our numerical 
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TABLE II: Comparison of hard- wall model to the e *W of Eq. ([3]) with Azq = 2.1; values in MeV 



(except for g piT1T ). 

Observable Experiment Hard wall \zq = 2.1 

m n 139.6±0.0004 [25] 139.6 139.6 

m p 775.5±0.4 [25] 775.3 777.5 

m ai 1230±40 [25] 1358 1343 

U 92.4±0.35 [25] 92.1 88.0 

/ p 1/2 346.2±1.4 [26] 329 325 

/a 1 / 2 433±13 [27, 28] 463 474 

3 p7r7r 6.03±0.07 [25] 4.48 4.63 



simulation indeed verifies this behavior. The F n (Q 2 ) prediction using the modified Xq(z) in 
Eq. ([2]) is never smaller than that of the original soft-wall model, and hence lies at least as 
far from the F n (Q 2 ) data. 

Second, we explore the consequences of modifying the background field as in Eq. 
repeated here for convenience: 

\ 2 Z 2 1 

e -*W = 6 0-1 (3) 

e A z o + e A 2 — 2 

As expected, the parameter A serves to interpolate between the hard-wall and soft-wall 
models. Specifically, for Xzq ~ 1 we find that the model prediction for F v is very similar to 
the original soft- wall result, while for Azo~2 the prediction is very close to the original hard- 
wall result: See Fig. [lj Beyond A^o~2 the global fit does not improve much; in particular, 
the F n (Q 2 ) fit is best for A^o = 2.1, but at the cost of f n = 88.0 MeV, while at A^o = 2.4 
the F 7T (Q 2 ) fit is not very good but one obtains / 7r = 92.7 MeV. That results mimicking the 
older models are obtained for such modest values of A^o is especially interesting because the 
hard- wall limit is formally recovered only when Xzq^oo, and the soft-wall limit is formally 
recovered only when Az^A^o, 1- 

As a comparison, we present results for both Xzq = 1 and Azo = 2.1 as well as the original 
soft- and hard-wall predictions in Tables [I] and [III respectively. Even though the result 
for A^o ~ 2 is very similar to original hard- wall result, its predictions for the positions of 
higher-order vector meson poles differ strongly from those of the original hard-wall model 
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Q (GeVj 

FIG. 1: Spacelike scaling behavior for F W (Q 2 ) as a function of Q 2 = —q 2 . The continuous line is the 
prediction of the original hard- wall model with 1 /zq = 323 MeV. The dotted line is the prediction 



of the original soft-wall model with K = m p /2. The crosses use e - *( z ) of Eq 
and the pluses use Xzo = l. The stars are from a data compilation from CERN 
from DESY, reanalyzed by Tadevosyan et al. 
the boxes 



3]) with Xzq = 2.1, 



the circ 



311 ] . the triangle is data from DESY 



es are 



3J, and 

3l| and diamonds [33j are from Jefferson Lab. Older data in the range 3-10 GeV 2 34] 



exist but have large uncertainties and are not presented here. 

(Table ITTTT) . The modified e - *^ results exhibit a closer match to a m 2 ~ n vector meson 
Regge trajectory, in contrast to the m 2 ~ n 2 scaling of the original hard-wall trajectory. 
From Table [TTT1 one finds in the hard-wall case that m n+1 — m n is stable around 1010 MeV, 
while in the modified background case, (m 2 +1 — m 2 ) 1//2 is stable around 1450 MeV. We 
conclude that the e~ x2 z2 tail behavior of Eq. ([3]) is already significant when A^ ~2, where 
the prediction for F V {Q 2 ) closely matches the hard- wall result. We have thus produced a 
model that works as well as the hard-wall model in fitting to the low-energy data, but that 
in addition reproduces the appropriate Regge trajectories expected from linear confinement. 



For completeness we also include the pion charge radius (r 2 ) results for both models. 
In the modified Xq model for B c = 0.96 and A c = 1.1, we find (r 2 ) = (0.485 fm) 2 . In the 
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TABLE III: Comparison of vector meson masses and decay constants (in MeV) for the hard-wall 
model and the model using e~*( 2 ) of Eq. ([3]) with \zq = 2.1. 



Original hard wall e *^ 2 - ) from Eq. ([3]) with \zq = 2.1 



771/2 rnl/2 

m p F p ' m p F p ' 

775.6 329 777.5 325 

1780.2 616 1608.1 528 

2790.8 864 2226.8 611 

3802.8 1089 2637.5 644 

4815.2 1300 2986.6 683 



modified e - * model for Xz = 1.0, (r 2 ) = (0.500 fm) 2 , while for Xz = 2.1 it is (0.576 fm) 2 . 
In comparison, the original soft-wall model gives (r 2 ) = (0.494 fm) 2 , the original hard-wall 
model gives (r 2 ) = (0.576 fm) 2 , and the experimental value is (r 2 ) = [0.672(8) fm] 2 [25]. 

Having addressed the two central issues of this paper, we now discuss an even more 
interesting conundrum, namely, the fact that all the fits described here are limited in how 
closely they approach the full set of F n (Q 2 ) and other low-energy data. All of them produce 
curves too shallow in Q 2 , leading in particular to a value of (r 2 ) smaller than experiment. We 
observed in Ref. that a better fit to F n (Q 2 ) could be obtained by loosening the fit to the 
other low-energy observables, particularly by letting /„. be somewhat smaller, a resolution 
that is not entirely satisfactory. However, as is well known 35], the functional behavior of 
F n (Q 2 ) for large Q 2 is determined by its partonic substructure: F n (Q 2 ) oc f 2 a s (Q 2 )/Q 2 . 
The AdS/QCD model used here contains only hadronic and no explicit partonic degrees of 
freedom, so it is perhaps not surprising that the F^i^Q 2 ) prediction has difficulty falling as 
steeply as the data for larger Q 2 . But are partons essential at the relatively small values 
(Q 2 < 5 GeV 2 ) presented in Fig. [Tp If not, then we are faced with the possibility that the 
pattern of chiral symmetry breaking developed in Ref. [5I and represented by the background 
field Xq(z) may require improvement. 
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IV. CONCLUSIONS 

In this paper we have shown, first, that the treatment in Ref. |24j of chiral symmetry 
breaking in the soft- wall model, although not formally satisfying all the boundary conditions 
of that model, nevertheless gives numerical results in agreement with those of models pos- 
sessing all the correct asymptotic behaviors. Indeed, we find that the soft-wall calculation 



of Ref. 24[ appears to represent a best-fit limit to such improved models. 

We have also developed a model that successfully interpolates between hard-wall and 
soft-wall models: It reproduces the desirable numerical features of the hard-wall model but 
also the Regge behavior predicted by linear confinement that motivates the soft-wall model. 
The fit to static observables as well as F n (Q 2 ) is as good as that from the hard- wall model, 
but the excited p masses now follow the desired pattern w? n ^n. 

Nevertheless, all of the models proposed here and in Ref. {24] share the properties that 
their predictions for F^{Q 2 ) are too shallow for larger Q 2 values, unless one relaxes the fit to 
one or more of the static observables (particularly /„-). We have argued that the most likely 
culprits for this discrepancy are either the application of a purely hadronic AdS/QCD model 
to a region (large Q 2 ) where partonic degrees of freedom become appreciable, or the need 
for an improvement to the pattern of chiral symmetry breaking introduced via background 
fields. Both directions are well worthy of further research. 



APPENDIX: NUMERICAL PROCEDURE 

Solutions to Eqs. (!6l)- (fT0|) are obtained numerically using subroutines provided by Nu- 
merical Recipes (NR) [36], and all terminology provided in this appendix is discussed in 
that well-known standard reference. These differential equations are solved using the shoot- 
ing method because the boundary conditions are defined on both ends of the integration 
range, generically labeled [21,22]. This procedure consists of several steps, which we now 
outline. Starting at the initial point 21 of the integration, suppose that N boundary condi- 
tions are to be specified but only n\ are provided. This leaves n2 = N— ri\ freely specifiable 
starting points, from which one can form an ^-dimensional vector V, the initial value of 
which is the first guess for the remaining boundary conditions. We then use the combined 
boundary conditions to integrate the differential equations and obtain 2/(22) (y being the 
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generic dependent variable label) at the endpoint £2 of the integration. Here one forms 
an ^-dimensional discrepancy vector F, which is the difference between the result y{x2) of 
the numerical integration and the specified boundary conditions at £2- The integration is 
performed via an adaptive stepsize Runge-Kutta (RK) method with constants provided by 
Cash and Karp. The solution is then obtained by optimizing the vector V such that the 
vector F converges to zero, by means of the globally-convergent Newton-Raphson method. 

As with any numerical procedure, this calculation has an inherent rounding error. There 
are two main sources for these errors, the RK integration and the Newton-Raphson method. 

1. The RK integration routine naturally provides in its result a 5th-order truncation er- 
ror (called jirA in NR) on the dependent variable. This truncation error scaled by 
some combination of the function y and its derivative dydx (using the terminology of 



Ref. 36J) evaluated at the initial point of each step in the integration determines the 
size h of the next step, but more importantly, drives the desired accuracy of our calcu- 
lation. In particular, the truncation error is obtained by taking the difference between 
the result calculated by adaptive stepsize RK (at 5th-order) and the result calculated 
by classical 4th-order RK. The small dimensionless parameter eps (fractional error) 
correlated to the truncation error is then defined by A = epsxyscal(i), where i is the 
index of the particular differential equation among those simultaneously solved that 
contributes the most to the truncation error (the "worst offender"). We choose the 
particular recipe for yscal in our calculation to be yscal(i) = \y(i)\ + \hxdydx(i)\ (see 
NR for other recipes). Lastly, we use eps = 10 -6 in our calculation. 

2. The Newton-Raphson method introduces several small parameters to assure conver- 
gence: 

• The small parameters ALF and TOLMIN are introduced to guarantee the proper 
convergence of the Newton-Raphson subroutine and therefore affect the final 
result of the calculation only indirectly. ALF functions as an indicator that 
the evaluated function has ceased to decrease appreciably, as determined by the 
Newton-Raphson method. In this case the subroutine backtracks along the New- 
ton direction until an acceptable rate of decrease is obtained. TOLMIN sets 
the criterion for identifying spurious convergence, in which case the subroutine 
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warns of the event and requires the introduction of a new initial point. These 
parameters are chosen to be ALF = 10 -4 and TOLMIN = 10~ 6 . 

TOLX sets the convergence criterion for 5x, which represents <5V in the shooting 
method algorithm. Its convergence criterion is set by TOLX = 10~ 7 . 

TOLF sets the convergence criterion on F values; TOLF = 10~ 4 . 



Further details can be found in Ref. 



36] 



All of these parameters are at least as small as (in most cases, appreciably smaller than) 
the measured uncertainties of the low-energy QCD observables to which we fit our models. 
As such, we can safely assume all the errors in our fits come primarily from experimental 
uncertainties that have been provided in the text, and we do not include further error 
analysis in our paper. Indeed, the chief figure of merit is the closeness of fit to the central 
values of certain hadronic quantities, as we now discuss. 

The meson masses m p , m^, and m ai are identified as poles in the solutions to Eqs. (jUJ)- 
(fID]l [see, e.g., Eq. ([16]) ]. One identifies pole masses by noting that the numerical solution 
near q 2 — > m 2 crosses a singularity and flips sign. Once a mass eigenvalue is identified, its 
individual bulk-to-boundary propagator is the solution to the relevant wave equation at that 
value of q 2 = m 2 , using the shooting method just described. With the normalization for ^5 
chosen, solutions to Eqs. (l6l)- (|T0l) at q 2 = are used to derive the decay constants f p , f^, 
and f ai [as in, e.g., Eq. (fl5l) ]. The overall quality of our fits is driven chiefly by (the central 
values of) the three best-measured meson quantities, m p , m n , and The other parameters 
in our models such as A c , B c , and A are then allowed to float freely. 
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